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Abstract 

Polaritons, in particular microcavity exciton-polaritons, have attracted much attention in 
recent years, as the phenomena of Bose-Einstein condensation and superfluidity have been ob¬ 
served for these quasiparticles. While the basic physics of these systems is well understood, 
there has been confusion over the connection of these systems to other phenomena, namely the 
Jaynes-Cummings Hamiltonian, Rabi flopping, and the vacuum Rabi splitting of atoms in a 
cavity. This paper reviews the basic theory of polaritons and shows these connections explicitly. 


1 Introduction 

Numerous recent works [I1I21I3I1 have reviewed the fascinating effects seen in experiments on 
Bose-Einstein condensation and superfluidity of microcavity exciton-polaritons. The field has 
progressed from simply observing evidence of Bose-Einstein condensation w to various effects 
such as Josephson junctions mi, phase locking of two condensates [9] , superfluid motion past 
a barrier m, and quantized circulation in a ring m- The lifetime of the polaritons in the 
cavities has steadily progressed, from about a picosecond in early experiments up to about 200 
ps in present systems m- 

Amidst all this activity, some basic confusions remain about the nomenclature for the system 
and the connection to other systems with light-matter coupling. For example, does the “Rabi 
splitting” of exciton-polariton systems have anything to do with Rabi oscillations or vacuum 
Rabi splitting? We will see here that indeed these are variations of the same phenomenon. 
Similarly, is the term “strong coupling” used in the polariton held in the same way as in other 
helds of optics? Again, we will see that this is the case. 

The basic concept of a polariton is simple. One starts with any oscillating dipole which can 
couple to the electromagnetic held. This can be, for example, an optical phonon in a solid, an 
exciton in a solid, or a two-level quantum oscillator consisting of two states in an atom. The 
coupling of the dipole oscillation to the electromagnetic held allows both radiation from the 
dipole or absorption of radiation by the dipole. 
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2 Phonon-polaritons 


The polariton effect is most easily introduced by examining the case of phonon-polaritons, when 
an optical phonon in a solid couples to the electromagnetic field. In this section we review the 
standard theory for phonon-polaritons. 

We start with the interaction energy of a dipole in an external electric field, i?dipoie = -qx-E. 
For a polarization field P = {N/V)qx, the total interaction energy is then 

Hint = -J d\ P-E. ( 1 ) 


The polarization, which is proportional to the local displacement, and the electric field can 
be written in terms of the phonon and photon operators. Dropping the vector notation for 
simplicity, the electric field is m 


E{r) = 



( 2 ) 


where Coo is the dielectric constant of the medium not counting the contribution of the optical 
phonons, and the polarization field is m 


^ y 2m.Ncoo 




c,P>^- + cle- 


ik-r'^ 


(3) 


where we have used creation and destruction operators cj, and Ck for the phonons, and we assume 
an optical phonon with constant frequency wq. This may be written more simply as 


P(r) = 


Ihe^ 

2V^ cl^o 




cte“- + 4e- 


ik-r'^ 


(4) 


where D = q'^N/eacixirV. Substituting these formulas for E and P into we obtain 

. hfl [uJk 1 
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The integral of the exponential factors over r gives us a 5k,w which eliminates one momentum 
sum, so that we have the total Hamiltonian 


( 6 ) 


H f -b nuQclck + ^ (ofcC.fc -b akcl - 4cfe - 4c4) 


This can be simplified by defining new operators which are linear superpositions of the 
creation and destruction operators which appear in the Hamiltonian. We define the new bosonic 
destruction operator ^k = ctkO-k + PkCk + IkoLk + dkc'_k its Hermitian conjugate for the 
creation operator, with the coefficients ak ■ ■ - Sk chosen such that 

k 
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which implies 
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The condition Q is equivalent to the matrix equation 
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( 8 ) 


Following the standard diagonalization procedure of setting the determinant to zero, we have 

u:^-oj^{col+u;l)+ojliol-n^ul=0. (9) 

This is equivalent to the standard phonon-polariton equation m 

( 10 ) 


e(oo) \uj^ ~ , 

with the bare photon frequency Wfc = ckj-^/e^ and wf. = w£ — using 


^ = e(oo)(wi - wf,). (11) 

and our definition of above. There are two positive-frequency solutions. 

When Wfc = wq, it is easy to show that in the limit wq 3> fl, the energies of the eigenstates 
are Tiuj = ?iwo ± hQ, and the corresponding eigenstates are ^k = {a-k T ick)j'/2- Away from the 
crossover region, when either Wfc ^ wq or Wfc S> wq, the eigenmodes correspond to nearly pure 
photon Ofc and phonon Cfc operators. 
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3 Exciton-polaritons: ensemble of two-level oscillators 
(Frenkel limit) 

We adopt the same picture as the phonon picture, but now imagine an ensemble of two-level 
electron oscillators at discrete locations i. The polarization of each oscillator is m 

^ (^{c\x\v)bl^bi^ + {c\x\v)*bl^b^^'j , ( 12 ) 

where we use b]^ and b^^^ for the fermionic creation and destruction operators of the electrons, 
and n = c,v represent the conduction and valence bands, respectively; since in the Frenkel 
picture, each oscillator has no spatial overlap with the others, we could equally well call these 
e and g for the excited and ground states of each oscillator, as in the case of an ensemble of 
atoms. 

We define the exciton operator in the Frenkel limit as 

= (13) 

^ i 

The inverse Fourier transform is 

= (14) 


so we have 


Pi = -i 


.qN 1 




V muJoVN ■ ^ 

where we write p^v = (c|p|u) = imLijQ{c\x\v) [T7]. We then have 

= - [ d^r P-E 


(15) 
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2eV V muJoVN ^ ^ 




h q^N 


/Wfe 


2 V emV V mhujQ y uiq 


(-PcvakC-k + PcyttkCl + PcvalCk - p*yalc[^^ . 


(16) 


As for the phonon-polaritons, we write = afeOj, -I- /dfcCj, -I- 7 fea_^ -I- 5kC_^. Then 


^ = E + ^^oClCk) + Hint (17) 

k k 
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and 


[^k,H] = huj^^ 
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(19) 


and F = pcv •\/2/m?iwo; |Fp is the oscillator strength. (We have used the approximation that 
the Ck operators are purely bosonic.) 

This is equivalent to the matrix equation 
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The determinant is 
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w4-w2(w2+w2)+w2w2-|ll|2w^ = 0. 

This is equivalent to the standard polariton equation 


2 c^hP (wfn — w2 

Wn — W^ 


( 20 ) 

( 21 ) 

( 22 ) 


with w|n = Wg — and w^ = ckj 

Cavity Polaritons. In a cavity we use w^ = (c/n)y^A:j| + fc^) ~ (c/n)(l + A:j|/2fc_L), with 

fc_L = tt/L. For the resonant case we set w^ = (c/n)7r/L = wg. Then the determinant equation 
is 


which has the solutions 


or 


- 2w^Wo + Wo(wo - lllp) = 0, 


(jj — loq dz u^Qri, 
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(23) 
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where the final approximation is valid when ^ wq- These correspond to the lower polaritons 
and upper polaritons so familiar in microcavity polariton optics. There will also be “dark,” 
uncoupled states at wq which lie in the middle of the upper-lower polariton splitting. This 
uncoupled Wq was in the case of the phonon-polaritons. 


4 Exciton-polaritons: Wannier picture 


The Wannier limit of excitons consists of the case in which the electron and hole (empty electron 
state) are no longer confined to the same oscillator, but instead, due to coupling between the 
oscillators, the electron and hole can migrate to different oscillators. In this case there will be a 
Coulomb attraction between the free electron and the hole, which effectively acts as a particle 
with positive charge. This leads to bound states of the free electron and hole that are exactly 
the same as the Rydberg bound states of a hydrogen atom, but with the energy scaled by the 
dielectric constant of the medium. 

For this calculation we use the interaction Hamiltonian in fc-space instead of in real space: 


Hi 


int 


= --E 

m 


k.k' 


2eVuJk 


{c\p\v) -f -f 


-{c\p\v)* (c 


afc4,/c'+fc^cfe' +44.fc'-fe^c,fc' 
The exciton creation operator in the Wannier case is [TS] 

k' 

The Fourier transform of the Is wave function is 






(26) 


(27) 


(j){k) = 


\/V (1 + aP'kP'Y ' 


(28) 


where a is the exciton Bohr radius. 

It is not easy to invert this to write the Hamiltonian in terms of the exciton operators. 
Therefore instead of an exact diagonalization, we write a matrix on the states |ex) = C^lO) and 
I phot) = a|,|0). The off-diagonal term is 


ex|H|phot) - (0|^(/)(fc/2 E l/p H 

k"> k" k' ^ ^ 

X 

{c\p\v) {aknbd+k''Kk' + ai,,bl,k'-k''K,k') 


+< 

(afe,.4,fc'+/c"^cfc' +4"4,/c'-fc"4fe') 

4lO) 


= -( 0 | ^ (/>(fc /2 - 

k'" 


(29) 

(30) 


(31) 
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m V 2eVuJk 
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(32) 


The sum over k' is 


^ ^ 2Trk'^dk'd{cos 9) 


(27r)3 

(27r)3 


2nk'^dk'd{cos9) 


2TTk'^dk' 


v^(l + a2|fc/2 + fc'|2)2 
SV 7ra3 

(1 + a^{\k‘^ + + kk' cos0))2 


(27r)3yo 2+|a4(fc2-4fc'2)2 + a2(fc2+4;./2) 


(33) 


Therefore 


(ex|i7|phot) = - 



(34) 


where 


n = 


emTra"’ 


tF 


and F is defined as above for the oscillator strength. 
We thus have the matrix 


(ex|7J|ex) 
(phot I i? I ex) 


(ex|77|phot) 
(phot I i? I phot) 


= n 


which has the determinant equation 


Wo ^o/‘ 


OJk 


(35) 


(36) 


WoWfe — uj{u}k + Wo) + W^ — - jnpWo/Wfe — 0. 


(37) 


This has the same behavior as the previous determinant equation (211 for the exact diagonal- 


ization, except near Wfc = 0 where it breaks down. See Figures 1 and 2. 


Note that the Rabi frequency (351 depends on the exciton Bohr radius a through the volume 


7ra3. The Frenkel limit can be viewed as the Wannier picture of excitons in the case when the 
exciton Bohr radius becomes equal to the unit cell size of the underlying crystal. This, in the 
Frenkel limit, the volume ira^ becomes the unit cell size a\ = V/N^ in which case the Rabi 


frequency (351 of the Wannier limit becomes exactly the same as (19) in the Frenkel limit. 
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Figure 1: Plot of the solutions of Eq. (37) for wq = 1, = 0.1. 


5 Comparison to atomic vacuum Rabi splitting 


As mentioned above, the Frenkel exciton limit is no different from the case of an ensemble of 
isolated atoms. We can therefore compare the case of an ensemble of atoms in vacuum to the 
above results with just a small change in notation. 

The standard Hamiltonian for vacuum Rabi splitting with a two-level atom is m 


where 


= V huJobl^bi^ + hujkala,, + V ng{akb\j3i^ + 


^ / ujqcP ^ ^ 1 j 

V 2‘heV V 2TieV 2^ emV ymTiwo/ 


(38) 


(39) 


and we have again used (p) = imuj{x). 
rescaling of the photon operator 


The Hamiltonian (381 is equivalent to (161 with the 




(40) 


assuming p^v is real, and dropping the two terms and a\b^J}\^, which is comparable to 

the approximation made in Section]^ which is valid when ojq ^ H. 

Agerwal [T5] found that the splitting is 


H = 2gVN, (41) 

which is exactly the same as found for our definition ([^ for the Rabi splitting of a two-level 
Frenkel exciton system. It really is the same system, namely an ensemble of independent two- 
level oscillators coupled only by the electromagnetic field, just solved in a different way. 
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It may at first seem as though the result (35) is not the same as the coupling (41), since the 
latter is proportional to '/N, where N is the number of atoms, while (35) is not proportional to 
the number of excitons. It is important to keep in mind, however, that the number of oscillators 
in an excitonic system is not given by the number of excitons, but by the number of atoms in 
the underlying medium. In the Frenkel case, the number of oscillators is exactly the number 
of atoms, while in the Wannier case, the effective number of oscillators is reduced by the ratio 
(ol/o)^, since each Wannier exciton is spread over many lattice sites. 


6 Final Remarks 


We have seen that the description of the interaction of electromagnetic field with all three types 
of electronic oscillation, namely Frenkel excitons, Wannier excitons, and exceptions of atomic 
states, can be described by the same formalism. The coupling term in each case can be viewed 
as the natural unit of frequency for the susceptibility of an ensemble of classical oscillators [50] 
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(42) 


where I4sc is the volume per oscillator. In each of the cases, “strong coupling” can be defined as 
the limit when the Rabi frequency 11 is much greater than 7 = 1/r, where r is the decay time. 

We have also seen that the standard approximation using a 2 x 2 matrix to represent the 
mixing of the photon states and electronic excitation is only valid in the limit when w ^ H. 
When n is comparable to w, sometimes called the “ultrastrong” coupling limit [5T|, the full 
Hamiltonian must be used. 

The use of the term “vacuum Rabi frequency” for the state splitting H is connected to the 
standard Rabi frequency. The standard Rabi frequency corresponds to the rate at which an 
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electronic oscillator will flip states in the presence of a classical driving field, and is given by 


— Eq, 


(43) 


q\pc 

mhuJo 

where Eq is the classical field amplitude. To get the vacuum Rabi frequency, we use the ampli¬ 
tude of the field which corresponds to the average electric field amplitude in a vacuum. From 
([^ we have 

huiQ 


which for a vacuum gives 


Then (43) becomes 






(44) 


(45) 


= 


q\p. 

mTiLOQ 



(46) 


which is the same as the coupling term we have used above. 

One of the important physical implications of these calculation is that the coupling energy for 
exciton-polaritons does not depend on the number of polaritons, to lowest order. The number 
of excitons is not analogous to the number of atoms in the atomic case; rather, the number of 
excitons corresponds to the fraction of atoms in excited states. At high density, the number 
of excitons can affect the splitting between the upper and lower polariton branches, giving an 
effectively reduced coupling fl, due to phase space filling, which reduces the oscillator strength, 
which comes about when we can no longer use the approximation that the exciton operators are 
bosonic, as we did above. This has been seen [53] to lead to a collapse of a microcavity system 
to weak coupling, leading to standard lasing. 
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